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ABSTRACT 

The star formation rate (SFR) of the Milky Way remains poorly known, with often-quoted values 
ranging from 1 to 10 M Q yr _1 . This situation persists despite the potential for the Milky Way to 
serve as the ultimate SFR calibrator for external galaxies. We show that various estimates for the 
Galactic SFR are consistent with one another once they have been normalized to the same initial mass 
function (IMF) and massive star models, converging to 1.9 ± 0.4 M yr _1 . However, standard SFR 
diagnostics are vulnerable to systematics founded in the use of indirect observational tracers sensitive 
only to high-mass stars. We find that absolute SFRs measured using resolved low/intermediate-mass 
stellar populations in Galactic H II regions are systematically higher by factors of ~2 — 3 as compared 
with calibrations for SFRs measured from mid-IR and radio emission. We discuss some potential 
explanations for this discrepancy and conclude that it could be allayed if (1) the power-law slope of 
the IMF for intermediate-mass (1.5 M Q < m < 5 M Q ) stars were steeper than the Salpeter slope, 
or 2) a correction factor was applied to the extragalactic 24 /im SFR calibrations to account for the 
duration of star formation in individual mid-IR-bright H II regions relative to the lifetimes of O stars. 
Finally, we present some approaches for testing if a Galactic SFR of ^2 M Q yr _1 is consistent with 
what we would measure if we could view the Milky Way as external observers. Using luminous radio 
supernova remnants and X-ray point sources, we find that the Milky Way deviates from expectations 
at the 1-3(7 level, hinting that perhaps the Galactic SFR is overestimated or extragalactic SFRs need 
to be revised upwards. 

Subject headings: Galaxy: fundamental parameters — galaxies: star formation — H II regions - 
ISM: supernova remnants 



1. INTRODUCTION 

Measuring the star formation rate (SFR) of a galaxy is 
like taking its pulse. The SFR quantifies the population 
of massive stars, which sets the composition and energet- 
ics of the interstellar medium (ISM) as well as the rate at 
which the ISM is dispersed through galactic winds and 
locked away into long-lived, low-mass stars. Our under- 
standing of galaxy activity and the growth of stellar mass 
throughout cosmic t ime depends critica l ly on measure- 
ment s of SFR (e.g., iMadau et all IT996: Bouwen s et al.1 
1201 If) . Galaxy evolution m odels require the SFR as a fun- 
damental input parameter (lLarsonll974HChiappini et al.l 
[19971: Kauffmann fc Haehneltil2000l:lde Rossi et al.H 2009). 
since the SFR sets both the level of massive star feedback 
and the rate at which gas and dust are depleted. 

The Milky Way is generally excluded from compara- 
tive studies of star-forming galaxies because we lack a 
unified, global picture of its structure and star forma- 
tion activity. The location of the Solar System inside 
the Galactic disk gives us a disadvantaged perspective 
when attempting to study the Milky Way as a whole 
galaxy. Sightlines through the disk suffer very high ex- 
tinction due to the dusty ISM, so SFR diagnostics which 
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depend upon optical/UV observational tracers (e.g., Ha) 
cannot be applied to the Milky Way, just as they can- 
not be applie d to external gala xies viewed nearly edge- 
on (see e.g., iKennicuttl Il998al and iCalzetti et all 120091 
for reviews of SFR indicators). Radio and far-infrared 
(far-IR) observations are unaffected by extinction and 
therefore h ave provided most of the Galactic SFR esti- 
mates (e.g.lMiljer fc Scalolll979t iGiisten fc Mezgerl ll982l; 
iMisiriotis et al.l 120061: iMurrav k, Rahman! 1201011 These 
wavelengths show great promise for providing conver- 
gence between extragalactic and Galac tic SFR mea- 
surem ents, because a calibration like the ICalzetti et al.l 
(20071 hereafter C+07) 24 /Ltm diagnostic for individual 
star-forming regions could be extended to the Galactic 
case. There remain, however, significant challenges stem- 
ming from uncertain distances to Galactic H II regions 
and confusion of multiple star-formation regions over- 
lapping along a given line of sight. 

In spite of the observational challenges created by ex- 
tinction, uncertain distances, and confusion, the Milky 
Way should be the ultimate SFR calibrator. The Galaxy 
offers the overwhelming advantage of high spatial resolu- 
tion. Infrared observatories have advanced to the point 
that it is now routine to resolve individual stars in young 
embedded clusters several kpc distant. This makes it 
possible to directly probe the young stellar populations in 
massive clusters ionizing Galactic H II regions, constrain- 
ing the stellar initial mass function (IMF) and stellar ages 
and investigating the detailed physical processes govern- 
ing star formation. The resulting detailed star forma- 
tion histories of Galactic H II regions and the associated 
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molecular clouds can then be calibrated against diffuse 
emission tracers, like integrated IR luminosity, that are 
used to study star formation in external galaxies. 

This paper seeks to answer three questions. (1) What 
are current estimates of the SFR of the Milky Way, and 
are they consistent with one another? Estimates have 
been presented in the liter ature for the past 35 years , 
ranging from ~1 M yr" 1 (jRobitaille fc WhitnevlfeoToh 
to -10 Mq yr" 1 (jGiisten fe Mezgerll 19821 ). Some of this 
variance is due to fundamental differences in the data, 
but, as we will demonstrate, systematics like the use of 
different IMFs or stellar models play a significant role. 
(2) How solidly are absolute SFRs known? To relate 
SFRs to other physical quantities like gas masses and 
stellar masses, it is not enough to calibrate SFRs in a 
relative sense; we need to know that a measured SFR of 
1 Mq yr _1 actually corresponds to 1 M yr _1 of material 
incorporated into stars. And finally, (3) Do Galactic SFR 
measurements give results consistent with extragalactic 
SFR diagnostics? Studies of the Galactic SFR generally 
use fundamentally different strategies, as compared with 
estimates of extragalactic SFRs; for example, Ha is not 
a useful diagnostic in the Galactic plane because of high 
extinction, and most external galaxies are too distant 
to study their resolved stellar populations. We would 
like to be certain that, if we measure 1 M Q yr _1 of star 
formation in an external galaxy, the same conditions in 
the Milky Way would also produce a measured SFR of 
1 M yr^ 1 . We identify two SFR diagnostics which can 
be applied in both the Milky Way and external galaxies: 
a galaxy's most luminous supernova remnant (SNR) and 
its population of X-ray point sources. 

Note that throughout this work we use the standard 
term "Galactic SFR" to refer to the SFR of the Milky 
Way, and "extragalactic SFRs" when referring to exter- 
nal galaxies. The subsequent sections of this paper are 
organized as follows: in §?? we define the IMF and stellar 
population synthesis models used as input assumptions 
in our SFR determinations, and in Sj3]we correct for these 
systematics in an effort to homogenize Galactic SFR es- 
timates as much as possible. In §J]we derive the SFR in a 
sample of Galactic H II regions from their resolved stellar 
populations and discuss the results in the context of sys- 
tematic errors inherent in absolute SFR determinations. 
In [J5]we present SFR tracers that can be self-consistently 
applied both to the Milky Way and to external galaxies. 
We summarize our conclusions in 

2. OUR ASSUMED INPUTS FOR CALCULATING SFRS 

As all SFR diagnostics only probe a subset of stel- 
lar masses, the IMF is a critical assumption in SFR 
estimates. The IMF is often expressed as </>(log m) = 
dN/d log m, the number of stars in the logarithmic mass 
interval between logm and logm + dlogm. We normal- 
ize our results to the commonly- used "b roken power-law" 
Kroupa IMF (iKroupa fc Weidnerl [20031 : used in the SFR 
calibrations of lKennicutt et al [120091 hereafter K+09): 

4>(\ogm) oc m~ r , (1) 

T = 0.3; 0.1 M Q < m < 0.5 M , (2) 

T = 1.3; 0.5 Mq <m< 100 Mq. (3) 

It is worth noting that the iChabrierl (|2005l) IMF yields 
total stellar masses that are practically identical to the 



Kroupa IMF when both functions are normalized to high 
stellar masses (see Figure [1]), implying that the SFRs cal- 
culated in this work should also represent SFRs assuming 
a Chabrier IMF. 

Most SFR diagnostics indirectly measure a specific 
electromagnetic signature of massive stars, usually Ly- 
man continuum (hydrogen-ionizing) photon flux, and 
then extrapolate this to a larger stellar population. 
Therefore, to accurately compare SFRs, we must use 
the same stellar population models and Lyman contin- 
uum photon rates in all cases. Here, we will use the 
stellar popul ation models employed in version 5.1 of 
Starburst9£0 ( Leitherer et alii 9991 : 1 Vazquez fc Leithererl 
I2005D and used in the SFR calibration of K+09, assum- 
ing solar metallicity. These m ake use of the recen t O sta r 
models of lSmith et all (120021) and [Martins et all (120051) 

According to iKennicutt et al.1 (j 19941 ) and iKennicuttl 
(1998a), a SFR of 1 M yr _1 produces a Lyman con- 
tinuum photon rate N c = 9.26 x 10 52 photon s _1 for the 
iSalpeterl ([T955) IMF (assuming a mass range of 0.1-100 
M ). The ratio of SFR to Lyman continuum photon rate 
is 1.44 times smaller for the Starburst99/Kroupa IMF 
formulation than for the Salpeter IMF (K+09), so we 
will adopt the following relationship between Lyman con- 
tinuum rate and SFR for the Starburst99/Kroupa IMF 
calibration: 

SFR [Mq yi- 1 } = 7.5 x 10~ 54 N c [phot s -1 ]. (4) 

This relationship will be used to compare Galactic SFR 
determinations throughout this paper. 

Equation [4] should be interpreted as the continuous 
SFR required to maintain a steady-state population of 
ionizing stars that produces an observed ionizing pho- 
ton rate N c . The "steady state" condition requires that 
the ionizing-star birth rate balance the death rate. The 
timescale tacitly assumed by Equation @] to compute 
SFRs is thus the lifetime of the ionizing stars. Sim- 
ple experiments with the Starburst99 code, assuming a 
constant SFR, demonstrate that on-going star forma- 
tion must last a minimum time of —8 Myr (the life- 
time of a late O star) to reach the steady-state value 
of iV c , and N c is within 5% of the steady- state value by 
~5 Myr, the lifetime of an early O star ()Bertelli et al.l 
ll994HMartins et alJl2005[ ). Unsurprisingly, O stars com- 
pletely dominate the ionizing photon flux from actively 
star-forming populations. We caution that the steady- 
state approximation underestimates SFRs for stellar pop- 
ulations with durations of star formation significantly 
less than O-star lifetimes. We will bear this issue in 
mind throughout §3.11 and return to it in §4.3.71 below. 

Where we determine SFRs from SN rates instead 
of Lyman c ontinuum photon ra tes, we assume Padova 
isochrones (|Bertelli et al.1 119941 ) and that stars in the 
mass range from 10-100 M explode as SNe. To de- 
termine SFRs in individual star-forming regions from 
intermediate-mass star counts, we also assume a Kroupa 
IMF and use stellar population ages from the literature, 
which gene rally are based o n pre-main-sequence (PMS) 
isochrones (jSiess et alJl200Q[ ). 

3. REVIEW OF MILKY WAY STAR FORMATION RATE 
ESTIMATES 

4 http:/ /www. stsci.edu/science/starburst99 
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3.1. Lyman Continuum Photon Rates 

ISmith et al.l (|1978[ ) established the "canonical" value 
for the Galactic SFR by calculating the total Lyman con- 
tinuum photon rate required to maintain the ionization of 
all Galactic giant H II regions then known from radio con- 
tinuum surveys. They reported N c = 4.7 x 10 52 phot s -1 
(shown by later studies to be a lower limit) based on the 
thermal free-free radio continuum flux emitted by the 
Galacti c H II regions i n thei r sample, using the formula- 
tion in iMezger et al.l (|1974D and correcting for a factor 
of ^2 in ternal absorption o f Lyman con tinuum photons 
by dust. ISmith et all (| 19781 ) then used a lSalpeterl (|1955f ) 
IMF (wit h limits of 0.1 and 100 M©), the massive star 
models of lPanag ia (1973), and an (unrealistically low) as- 
sumed average duration of star formation in a Galactic 
radio H II region of 5 x 10 5 yr to transform from Lyman 
continuum photon rate to a Galactic SFR of ^5 M© yr -1 
(factor of ^3 reported uncertainty). Their measurement 
of jV c translates to a Galactic SFR of 0.35 M© yr -1 with 
the Starburst99/Kroupa IMF calibration used here (Ta- 
ble [1]) . This strikingly low value can largely be explained 
by the incompleteness of their H II region sample, but is 
also contributed to by the inappropriatcncss of the Star- 
burst99 steady-state calibration (Equation |4j as applied 
to the young H II regions in the Smith et al. sample (this 
will be discussed further in N4.3[) . 

Building upon the work of ISmith et ail (|1978l ). 
IGiisten fe Mezgerl (|1982l ) made a new estimate of the 
Galactic SFR from radio H II regions, taking into ac- 
count the thermal radio emission from "extended low- 
density" H II regions and including Lyman continuum 
photons escaping from giant H II regions. They mea- 
sured N c — (2.7 ± 0.8) x 10 53 phot s _1 , nearly six times 
the Smith et al. value, and reported a Galactic SFR of 
11 M Q yr -1 , which is clearly not six times the S FR es- 
timate of 5 M yr -1 made bv ISmith et al.1 (| 19 781 ). This 
discrepanc y is explained by the adoption of the lognor- 
mal IMF of iMiller fe Scald (| 19791) . different limits on the 
IMF mass range (from 0.1 to 60 M ff ), and different st el- 
lar models (see Appendix A of lGiisten fe Mezger]|1982l for 
a description of their conversion from Lyman continuum 
photon rate to stellar mass). In our homogenized Star- 
burst99/Kroupa IMF calibration, the IGiisten fe Mezgerl 
( 1982) measurement of N c translates to a SFR of 2.0±0.6 
M© yr -1 (Table [1]). The incorporation of low-density 
H II regions means that this measurement included the 
contribution to i onization provided by older Galactic O 
stars omitted by ISmith et al.l (|1978[ ) , hence the steady- 
st ate calibration should be appropriate here. 

IMez ger (1987) calculated a slightly lower value for the 
Lyman continuum photon rate of N c = (2.1 ± 0.6) x 
10 53 phot s , discrepant from their earlier value due 
to their use of R© = 8.5 kpc for the distance to the 
Galactic center (as opp osed to R© = 10 kpc used by 
IGiisten fc Me^gen [198211 . This discrepancy underscores 
the point that Galactic SFR measurements are often de- 
pendent on models of Milky Way structure and star for- 
mation therein. Here, we are not accounting for differ- 
ences in the assumed Milky Way model. They also used a 
modified Millcr-Scalo IMF with a less steep slope for the 
highest mass stars, and find a SFR of 5.7 ± 1.6 M© yr . 
Using our Starburst99/Kroupa IMF calibration, we find 
a corresponding SFR of 1.6 ± 0.5 M© yr -1 . 



iBennett et al.l (|1994l ) used Cosmic Background Ex- 
plorer ( COBE) observations of N II 205 /im emission to 
measure the rate of Lyman continuum photons in the 
Milky Way, and correcting for dust absorption and pho- 
ton escape, th ey found N n = (3.5 ± 1-8) x 10 53 phot s _1 . 
Subsequently, iMcKee fe Williamsl ((1997) used the same 
COBE data and slightly different assumptions to mea- 
sure a Lyman continuum photon rate of N c — (2.6 ± 
1.3) x 10 53 phot s _1 . These two estimates of N c agree 
with one another to within the 50% uncertainties quoted 
by the authors, and are also in excellent agre ement with 
the radio data from IGiisten fe Mezgerl (|1982|) . 

In a r ecent update to the origin al Galactic SFR es- 
timates, iMurrav fe Rahman! (|2010( ) used Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) observations of free- 
free emission to measure a dust-corrected Lyman contin- 
uum photon rate of N c = 3.2 x 10 53 phot s — 1 . This is in 
good agreement with previous measurements of iV c and 
corresponds to a SFR of 2.4 M© yr" 1 (Table [J). The 
COBE and WMAP studies utilize the integrated light of 
the Milky Way and average over its entire stellar popula- 
tion; therefore, the steady-state Starburst99 calibration 
is also appropriate to these studies. 

3.2. Massive Star Counts and Supernova Rates 

Measurements of the core-collapse SN rate in the Milky 
Way have also been used to constrain the SFR. A SN 
rate of 0.85 events per century is expected for a SFR of 1 
M© yr" 1 , assuming a Kroupa IMF and a SN progenitor 
mass range of 10-100 M© (calculated using Starburst99 
models with t he parameters discussed in §2). 

iReedl (|2005| ) surveyed the B2-03 type stars (>10 M ) 
in the local Solar neighborhood (within 1.5 kpc) and as- 
sumed a model of the Milky Way to extrapolate these 
massive star counts to the entire Galaxy. He estimated a 
steady-state SN rate of 1-2 events per century, averaged 
over the lifetimes of massive stars (^20 Myr for 10 M©). 
In our Starburst99 /Kroupa IMF calibration, this trans- 
lates to a SFR of 1.2-2.4 M© yr" 1 . 

iDiehl et al.1 ([2006) used gamma-ray measurements of 
radioactive 26 Al to constrain the core-collapse SN rate 
in the Milky Way to 1.9 ± 1.1 events century" 1 , av- 
eraged over a timescale set by the half l ife of 26 Al 
(~ 7.2 x 10 5 yr). They assumed a IScaTol (fl998l) IMF 
with r = 1.7 and a SN progenitor mass range of 10-120 
Mp) . Assuming massive s tar yields of 26 Al as modeled by 
iLimongi fe Chieffil ()2006l ). the yields for a Kroupa IMF 
are a factor of 1.12 higher than for the Scalo IMF. This 
means that for our normalized assumptions, the SN rate 
from radioactive 26 Al is 1.7 ± 1.0 events century -1 , con- 
sistent with the results from Reed (2005). This translates 
to a SFR of 2 ±1.2 M© yr" 1 . 

3.3. Infrared Diffuse Emission and Point Sources 

iMisiriotis et a l. (2006) constructed a model of the spa- 
tial distribution of stars, gas, and dust in the Galaxy 
constrained by COBE observations of the Galactic IR 
emission. They measured the 100 /zm luminosity of 
the Milky Way and, u sing the conversion described in 
IMisiriotis efahl (|2004D . calculated a SFR of 2.7 M© 
yr- 1 . However, they used a Salpeter IMF (0.1-100 M ) 
and lFioc fc Rocca-Volmerangd (Il997) stellar population 
models. We attempt to correct for the difference in IMF 
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TABLE 1 

Normalized SFR Estimates for the Milky Way 



Method 


Observation 


SFR 


Reference 






(M© yr" 1 ) 





Ionization Rate 
Ionization Rate 
Ionization Rate 
Ionization Rate 
Ionization Rate 
SN Rate 
SN Rate 
Dust Heating 
Ionization Rate 
YSO Star Counts 
YSO Star Counts 



Radio Free- Free 
Radio Free- Free 
Radio Free-Free 
COBE N II 205 fim 
COBE N II 205 fj,m 
O/B Star Counts 
INTEGRAL Gamma-ray 26 Al 
COBE FIR continuum 

WMAP Free-Free 
Spitzer IR Photometry 
MSX IR Photometry 



0.35 Smith et al. (1978) 

2.0 ±0.6 Giisten fc Mezge77i982) 

1.6 ±0.5 Mezger (1987) 

2.6 ± 1.3 Bennett ct al. (1994) 

2.0 ±1.0 McKee fc Williams7l997) 

1.8 ±0.6 Reed (200*57 

2.0 ±1.2 Diehl et al. (2006) 

1.9 ± 0.8 a Misiriotis et al. (2006) 
2.4±1.2 a Murray fc Rahman7l5l0) 

1.1 ±0.4 Robitaille fc Whitney (2010) 
1.8 ±0.3 Davies et al. (2011) " 



Assumed 50% uncertainty, as no error is quoted by the original author. 



by dividing the SFR by 1.44 (as above, see Equation @| , 
assuming that the Lyman continuum photons and the 
non-ionizing UV photons which produce FIR emission 
scale together (a good assumption, as the disparity be- 
tween the Salpeter and Kroupa IMFs occurs at low stellar 
masses). Typically, galaxy- wide mid-IR measurements of 
SFR are sensitive to star formation over longer timescales 
than other SFR diagnostics (^100 Myr; C+07), as even 
non-ioni zing B stars can contribute significantly to dust 
heating (|Cox et all 119861 ). We note that FIR diagnos- 
tics of the SFR include a host of critical assumptions, 
like what fraction of stellar light is absorbed by dust 
and how much older stellar populations contribute to 
dust heating. Still, th e estimate of Galactic SFR from 
iMisiriotis et al.l (|2006[ ) is consistent with the others de- 
s cribed here, at 1.9 M w yr" 1 ( Table P . 

iRobitaille fc Whitnevl (|2010l ) produced Galactic popu- 
lation synthesis models of young stellar objects (YSOs) 
that still retain dusty circumstellar disks and/or natal 
envelopes, assuming Kroupa IMF, a simple prescription 
for accretion, and a YSO lifetime of 2 Myr. These 
pop ulation synthesis mode ls were then used to simulate 
the IRobitaille etaLl (|2008D catalog of YSOs identified 
via mid-IR excess emission from the Spitzer Space Tele- 
scope GLIMPSE (Gala ctic Legacy Infrared M id-Plane 
Survey Extraordinaire; iChurchwell et al.l 120091 ) surveys 
of the Galactic plane. From comparing the simulations 
to the catalog, IRobitaille fc Whitnevl (|2010l) derived a 
SFR of 0.68-1.45 M yr" 1 . IRobitaille fc Whitnevl ff2T)10h 
note that their quoted range of SFRs includes only 
the uncertainties on t he total number of YSOs in the 
IRobitaille et al.l (j0 08) catalog. 

iDavies et al.l ()201lD modeled massive (m > 8 M©) 
YSOs found in the Red MSX Sources (RMS) Survey, 
using the data to find the best fits to basic YSO pa- 
rameters such as bolometric luminosity and mass. This 
study assumed a massive YSO lifetime of ^0.1 Myr and 
a Kroupa IMF, resulting in a Galactic SFR of 1.5-2.0 
M© yr" 1 . 

3.4. Galactic SFR Determinations Compared 

We see remarkably good agreement in Tabled] between 
estima tes of the Galactic SFR (with the exception of the 
oldest, iSmith et all 119781 which was clea rly incomplete 
and superseded bv lGiisten fc Mezgerf l982), after we have 
taken care to normalize assumptions about the parent 



stellar populations. The estimates clus ter around an av- 
erag e of 1.9 ± 0.4 M w vr" 1 (excluding ISmith et alii 19781 
and iGiisten fc Me zger 1982, which were both improved 



upon bv lMezgerlll987). The most s i gnific ant outlier is the 
result from IRobitaille fc Whitnevl (|2010D . We note that 
their Galactic SFR could be regarded as a lower limit, 
due both to issues of completeness and to the 2 Myr 
lifetime of the YSO phase used in their preliminary pop- 
ulation synthesis models. This age limit is likely an over- 
estimate, because it is based on circu mstellar disk life- 
times for low-mass PMS stars, while the Robitaille et all 
(2008) catalog is sensitive primarily to intermediate- 
mass YSOs (masses >3 M©). There is mounting evi- 
dence that the optically thick inner disks of intermediate- 
mass Y SOs are rapidly destroyed over timescales 
<1 Myr (IPoyich fc WhiTn^[2TJTfl IKoenig fc Allenll2011l: 
IPovich et al.H201lD. Such a c orrect ion in a future update 
of the Robitail le fc Whitnevl (|2010D population synthesis 
models could easily provide the factor of ~2 increase nec- 
essary to bring the Galactic SFR from intermediate-mass 
YSO star counts into agreement with the other estimates. 

We stress that our homogenized estimate is signifi- 
cantly different from the Galactic SFR of 4-5 M© yr -1 of- 
ten quoted in the literature, and this disparity is due only 
to differences in the IMF and stellar models used. Inves- 
tigators hoping to compare properties of the Milky Way 
with oth er galaxies (e.g., the lumin osity functions of star 
clusters; lHanson fc Popescul 120081 ) must first calibrate 
the SFRs to the same scale. The Starburst99/Kroupa 
IMF calibration is likely to be used often in the future, in 
which case the corresponding Galactic SFR is 2 M© yr -1 . 

4. CALIBRATING SFR DIAGNOSTICS AGAINST DIRECT 
STAR COUNTS 

We have demonstrated that various measurements of 
the Galactic SFR using myriad observational tracers con- 
verge to a consistent value, provided that systematic dif- 
ferences are minimized by adopting a uniform IMF and 
stellar models. We have made no assertions that our 
assumed inputs for calculating SFRs (§??) are the cor- 
rect ones; they were chosen simply to match the Star- 
burst99 stellar population synthesis models widely ap- 
plied to SFR determinations in external galaxies. If these 
assumptions were incorrect, then relative SFR compar- 
isons would be unaffected, but potentially large system- 
atic errors would be introduced into absolute SFR deter- 
minations. 
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Fig. 1. — Simulated Galactic-field IMF (solid curve; 

IKroupa & Wcidncr 2003) a nd l ognormal+power-law IMF 

(dash-dotted curve; IChabrierl [2005) normalized to the same 
relative numbers of early O-type stars (m > 40 M©). Dotted 
curves show an error envelope on the Chabrier IMF based on 
an uncertainty of 0.3 in the value of the power-law slope for 
m > 1 Mq. Five mass ranges (numbered) relevant to multi- 
wavelength determinations of the IMF and SFRs are indicated 
schematically by brackets and discussed in the text. 

While the gene ral broken-p ower-law (or 
lognormal+power-law; IChabrierl 120031) form of the 
observed IMF is remarkably consistent across stellar 
populations in widely disparate environments, variations 
have been reported both in the power-law slopes and 
in the l ocations of the break points wher e the slopes 
change (|Kroupal 120011: iBastian et all 120101) . Although 
these variations have low significance given the typically 
large uncertainties inherent in IMF measurements 
(|Bastian et all I2010j ). even an insignificant deviation 
from the default Salpeter-Kroupa slope can have sig- 
nificant consequences. The center-of-mass of the IMF, 
defined as the point where the integrated stellar mass 
is divided equall y between hig her and lower masses, 
falls near 1 M flKroupal 120011 ). If we estimate SFRs 
by extrapolating over an IMF that is pinned only 
to massive stars, we maximize the impact of IMF 
variations because such massive stars arc farthest from 
the center-of-mass. This effect is illustrated in Figure Q] 
in which th e logn ormal+power-law (T = 1.35) IMF 
of IChabrierl ((2005) is plotted alongside the standard 
Kroupa IMF. Both IMF models have been normalized 
to the same relative numbers of very high-mass stars 
(early O stars with m > 40 Mq, which dominate the 
bolometric luminosity and Lyman continuum photon 
flux in observable extragalactic star- forming regions; 
iMartins et al.l 120051 ). With this normalization, the 
Chabrier IMF yields 15% fewer stars (defined henceforth 
as A^pop) but 3% greater integrated stellar mass (Af pop ) 
above the hydrogen-burning limit (m > 0.08 Mq) com- 
pared to the Kroupa IMF. The insignificant difference 
in Mp p means that the choice between these IMFs is 
unimportant for deriving SFRs from unresolved stellar 
populations in H II region s. However, an uncertainty of 
0.3 on the value of T (e.g. iMuench et all 120021; IChabrierl 
2005) increases iV pop by a factor of 2.3 and M pop by a 
factor of 1.6 on the high side envelope (upper dotted 
curve in Figure [J) . 



High-resolution, wide-field, multi-wavelength imaging 
observations enable us to derive SFRs in bright Galac- 
tic H II regions (the nearby analogs of observable extra- 
galactic star-forming regions) using direct detections of 
the resolved low- and intermediate-mass stellar popula- 
tions. Such an approach samples the center of the IMF, 
greatly reducing the uncertainty of extrapolation (Fig- 
ure [I} while avoiding altogether any dependence on mod- 
els of massive star properties. Hence the following anal- 
ysis provides an independent estimate of the systematic 
errors on absolute SFR calibrations. 

4.1. The Galactic H II Region Mil as a Case Study 

4.1.1. The X-ray Luminosity Function: SFR from Star 

Counts 

W e begin by considering the well-studied, nearby (2.0 
kpc; IXu et aljlMTlT ) Galactic H II region M17. Because 
of its relatively high obscuration {Ay > 8 mag), bright 
nebulosity, and overwhelming field-star contamination 
produced by its location in the inner Galactic plane 
([/, b] = [15?1, — 0?7]), establishing a reliable membership 
catalog for NGC 6618, the massive cluster ionizing M17, 
has proven prohibitively difficult using optical/IR p ho- 
tometry alone (see discussions in iBroos et al.l [20071 and 
iPovich e"t~a l. 2009). These problems have been overcome 
by employing high-resolution imaging spectroscopy from 
the Advanced CCD Imaging Spectrometer (ACIS) of 
the C handra X-ray Obse rvatory to identify cluster mem- 
bers (|Broos et al.f l2007f ). Low- and intermediate-mass 
PMS stars exhibit powerful convectively-driven magnetic 
reconnection flaring activity, hence they are 10 2 -10 3 
times more luminous in X-rays than older Galactic field 
stars, and their X-ray luminosity Lx is broadly corre- 
lated with bolometric lumin osity and hence stellar mass 
(jPreibisch fc Feigelsonl[2005T ) . Hard (>2 keV) X -ray pho- 
tons are relatively unaffected by interstellar absorption, 
and while optical/IR photometry in bright H II regions is 
plagued by diffuse nebular emission, ACIS point-source 
detections generally do not suffer si gnificant reduction s 
in sensitivity due to diffuse emission flBroos et al.l l2007). 

The stellar X-ray luminosity function (XLF; the distri- 
bution of stars as a function of Lx) itself provides a good 
statistical tracer of the underlying IMF i n young clusters 
(|Getman et alj l2006t IWang et all 120071) . The approxi- 
mate range of masses sampled by observed XLFs of vari- 
ous massive clusters is indicated in Figured] In addition 
to PMS stars, X-ray observations are highly sensitive to 
massive stars, which produce strong X-ray emission via 
microshocks in their stellar winds and a variety of more 
exotic mechanisms ([Gagne et al.l 1201 ID . Intermediate- 
mass A- and B-type stars on or near the zero-age main se- 
quence lack both magnetic dynamos and po werful stellar 
wind s, hence they are generally X-ray quiet (jPovich et al.1 
120111 and references therein). 

To estimate -/V pop from t he XLF of a given c luster, 
the simplest approach (e.g., iFeigelson et al.l [20 111 ) is to 
use as a calibrator the XLF of the ONC, which contains 
839 lightly-obscured (X-ray median ene rgy < 2 keV), 
low-mass members (jGetman et al.l f2005t ) . The XLF in 
the target cluster is compared to the calibrator ONC 
XLF to derive a scale factor, which is then multiplied 
by iVp Op (0NC), th e population of the ONC abov e the 
brown-dwarf limit (Hillenbr and fc Hartmannlll998l) . 
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Estimates of -/V pop from XLF scaling assume several 
fundamental similarities between the target cluster and 
the ONC: (1) the form of the underlying IMF, (2) the ra- 
tio of obscured/unobscured cluster members, (3) the cir- 
cumstcllar disk fraction of the cluster membership (Lx 
is lower for protostars and YSOs with disks compared to 
diskle ss PMS stars of the same mass; see iPovich et al.l 
120111 and references therein) , and (4) the transformation 
between XLF and IMF. With the possible exception of 

(1) , none of these assumptions is generally valid. As- 
sumptions (3) and (4) could hold if the age of the tar- 
get cluster is sufficiently close to that of the ONC, but 
in practice ages are difficult to determine, and in any 
case it is not yet known how much the target cluster 
could differ from the ONC age before producing a signif- 
icant impact on the XLF scaling. To avoid assumption 

(2) , a "two-component" approach can be employed, in 
which the cluster XLF is divided into lightly-obscured 
and heavily-obscured sub-populations, and the calibra- 
tor ONC XLF is scaled independently to the XLF of 
each sub-populat ion (this second approach was used by 
iBroos et al.ll2007l ). 

In M17, the XLF was measur ed from 886 source s 
detected by Cha ndra/ ACTS (|Broos et al.l I2007D . 
IPovich et all (J2009D showed that Ap 0p as reported 
bv IBroos et al l (|2007[ ) for NGC 6618 translates to a SFR 
based on s tar counts of SFRsc = 0.008-0.01 M yr" 1 , 
using the iMuench et~aT1 (|2002j) IMF to convert A pop 
to Mp 0p . Adopting instead our standard Kroupa IMF 
an d an upper envelo pe of 1.5 Myr (see discussion 
in IPovich et al.l 120091 ) for the age spread of the X- 
ray-detected stellar population in NGC 6618, SFRsc 
becomes 0.004-0.006 M Q yr- 1 . 

In the following two subsections, we estimate the 
M17 SFR using independent tracers of the massive stel- 
lar population and compare the results to SFRgc ~ 
0.005 M Q yr- 1 . 

4.1.2. The Thermal Radio Continuum SFR Diagnostic in 

Mil 

Historically there has been some disagreement between 
different measurements of the M17 radio continuum flux 
densi ty. Adopting ffschz = 784 Jy and SlsGhz = 609 Jy 
from IPovich et al.l (|2007| ) and using the same conversion 
from thermal continuum (free-free em i ssion) to Lyman 
continuum photon rate as iSmith et al.l (|1978f ) , we calcu- 
late N' c = 3.07 x 10 50 s" 1 at 5 GHz and 2.66 x 10 50 s" 1 
at 15 GHz. The difference reflects a ~15% uncertainty 
on N' c . 

Summing over the known O sta r s in N GC 6618 and 
using the models of iMartins et al.l j2005) to obtain the 
Lyman continuum photon emission rates as a function 
of spectral type, the combined Lyman continuum pho- 
ton rate emitted is N c = 3 x 10 50 s" 1 . This is m sur- 
prisingly good agreement with N' c derived above, given 
that some fraction of the Lyman continuum photons are 
expected to be absorbed by dust before contributing to 
the ionization of the H II region, and some fraction of 
the Lyman con tinuum flux likely e scapes from the neb- 
ula altogether (|Povich et al.l 120071 ) . Our calculation of 
N c largely neglected the contribution of stars later than 
08 to the ionization, but the inclusion of 08-B3 stars 
would only increase N c by a few percent. There may 
yet be new major ionizing stars waiting to be discov- 



ered in M17 ()Povich et al.l 120081 reported the discovery 
of a new candidate 05 V star), but it is unlikely that 
a significant fraction of the ionizing stellar population 
in M17 has been missed. We have accounted for known 
massive binaries among the majority of the principal ion- 
izing stars (types 6 and earlier) in our calculation of N c 
(jPovich et al.l 12009;). We note that the total bolometric 
luminosity of the known O stars is Lob ~ 5.8 x 10 6 L , 
in very good agreement with the luminosity in the global 
spectral energy dis tribution (SED) of the M17 nebula 
(jPovich et al.ll2007D . 

Absorption of Lyman continuum photons by dust has a 
negligible impact upon the ionization of M17. Somewhat 
more important is the escape of Lyman continuum pho- 
tons into th e diffuse ISM thro ugh the blow-out in the 
H II region ([Povich et al.l I2007D . I n the most probabl e 
nebular geometries, N' c /N c > 0.9 (jPovich et all 120071) . 
Given the uncertainty in N' c , it is not unreasonable to 
assume N c = N' c = 3 x 10 50 s _1 for this highly em- 
bedded H II region. Using our normalized conversion in 
Equation El this translates to SFR// = 0.0022 M Q yr -1 
traced by free-free emission, which is lower than SFRgc 
by a factor of >2. In £ 14.31 we will discuss some possible 
explanations for this discrepancy. 

4.1.3. The 24 [im Mid-IR SFR Diagnostic in Mil 

According to C+07, their Spitzer 24 fim SFR calibra- 
tion "is appropriate for metal-rich H II regions or star- 
bursts" where the radiation field is dominated by young 
massive stars. M17 is a highly embedded Galactic H II 
regi on that radiates > 90% of its luminosity in the mid- 
IR ()Povich et al.1120071) , hence it meets these criteria. We 
can therefore use Ml 7 to test the applicability of the 
C+07 extragalactic mid-IR SFR diagnostic to a nearby, 
resolved H II region. 

M17 was observed using the Multiband Imaging Pho- 
tometer (MIP S) on board Spitz er as part of the MIPS- 
GAL survey (|Carev et all I2009D . but because the H II 
region is bright enough to completely saturate even the 
short MIPSGAL exposures at 24 /im, there is no direct 
measurement of the luminosity in the MIPS 24 fim band, 
L 2 4, that can be used as an input to the C+07 SFR re- 
lation (their Equation 9) Instead, we interp olate the 
global SED presented by IPovich et all (|2007l) over the 
range of 20 to 60 /im and compute the 24 fim flux den- 
sity, F„(24) (see Figure Eh). This gives 

L 2i = vL v (2A) = 4nd?vF v (2A) = 6.6 x 10 39 erg s _1 , 

(<25% uncertainty, dominated by the uncertain dis- 
tance d) and, following C+07, SFR 24 =SFR// = 0.0022 
Mq yr -1 , just as we derived from radio free-free emission 
above. We therefore find that SFRs based on indepen- 
dent indirect tracers of Lyman continuum photon rate 
are in agreement, but both are lower than SFRgc by fac- 
tors of >2. We note that the C+07 calibration uses the 
Kroupa IMF and Starburst99 models as we have used 
throughout this paper, including in our calculation of 
SFR//. As cautioned in §??, M17 may be too young for 
the steady-state approximation to hold. We will discuss 
the consequences of ignoring this caveat in ^4.3. 71 

The C+07 L24 SFR diagnostic is a statistical relation 
based on 179 observed H II region "knots" in 23 differ- 
ent high-metallicity galaxies, and its disagreement with 
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TABLE 2 

SFR Data for 8 Galactic H ii Regions 





Alternate 


d 


Earliest 




T~SF 


log L 2 A 


logA^ 


SFR SC 




Name 


Name 


(kpc) 


known star(s) 


(M ) 


(Myr) a 


(erg s — 1 ) b 


(phot s" 1 ) 


(10- 3 M yr" 1 ) 


Refs. d 


Carina 


NGC 3372 


2.3 ±0.1 


LBV,WNL,02 I 


43000 ± 5000 


5 


40.04 ± 0.04 


50.9 ± 0.1 


8.6 


1,2 


RCW 49 e 


Wd 2 


4.2 ± 1.5 


4xWNL,03 V 


>4000 


3 


39.98 ±0.35 


50.5 ± 0.3 


>1.3 


3,4 


M17 


NGC 6618 


2.0 ±0.1 


4x04 V 


8000 ± 1000 


1.5 


39.82 ±0.08 


50.5 ± 0.1 


5.3 


5,6,7 


NGC 6357 


Pismis 24 


1.7 ±0.5 


03.5 I, 03.5 III 


6500 ± 1300 


2 


39.24 ±0.20 


49.6 ± 0.3 


3.2 


8,9,10 


RCW 38 




1.7 ±0.2 


05.5 V 


1200 ± 300 


1 


39.03 ±0.10 


49.8 ± 0.1 


1.2 


11 


W3 Main 




2.0 ±0.07 


Several ~05 


3250 ± 1300 


1.5 


38.93 ±0.03 


49.45 ± 0.05 


2.2 


12 


Rosette 


NGC 2244 


1.3 ±0.1 


04 V 


1300 ± 200 


2 


38.63 ±0.10 


49.8 ± 0.1 


0.65 


13,14 


ONC 


M42 


0.41 ±0.04 


07 V 


1600 ± 100 


3 


38.48 ± 0.08 


48.77 ±0.09 


0.53 


15,16 



a The star formation timcscale tsf is defined as the age of the oldest stars included in the IMF scaling that yielded A/ pop . Hence SFRsc is the 
time-averaged SFR over tsf, and the actual SFR during any short bursts could be significantly higher, 
k L24 = ALa is the equivalent mid-IR continuum luminosity in the MIPS 24 /im bandpass. 

c N' c is the Lyman continuum photon rate required to maintain the ionization of the radio H II region, following the notation and method of calculation 
employed bv lSmith et al.l U978I) . 

d 1 HSmith&rBr ooks (2007), 2 = IPovich etTaE! (120111), 3 = ITsuihnoto et alj J2007I) , 4 HAscenso et al.l (t2007T) . 5 =[Povich et al.l J2007ft . 6 =IBroos et aTl 
Povich et al. (2009), 8 = Wang et al. (2007), 9 = Russcil et al. (2010), 10 = Maiz Apcllaniz et al. (2007), 11 = Wolk et al. (2006), 12 = 
ICelnikl U985D , 14 =IWang et all fl2008fiT 15 = [Hillenbrand k. Hartmannl fl 19981) , 16 = IFeigelson et al.l J2002IK 



J2007I). 7 

Fcigelson & Townslcy (2008), 13 



and SFRsc ar e lower limits for RCW 49 because the 2.8 kpc distance dAsccnso et al. 2007) assumed for the IMF scaling is at the near end of 



the wide range of distance estimates reported in the literature. 

SFRsc in a single, Galactic H II region is hardly a cause 
for alarm. Perhaps M17 is simply an outlier from the 
trend, or perhaps there is some fundamental difference 
in measureme nt between the Gal actic and extragalactic 
observations. iPovich et al.l (|2007T > measured the mid-IR 
luminosity of M17 using a large aperture chosen to match 
the linear extent of the H II region and photodissociation 
region, while the H II knots in distant galaxies were gen- 
erally smaller than the apertures used by C+07 to extract 
their fluxes. Hence, there could be an offset between the 
mid-IR luminosity measured for Ml 7 and the luminosi- 
ties of the extragalactic H II knots due to the drastically 
different filling factors and background subtractions in- 
volved. The impact of such potential biases can be eval- 
uated by hypothetically locating M17 among the C+07 
extragalactic sample (their Figure 9). For each H II knot 
in their sample, C+07 measured "luminosity surface den- 
sities" (LSDs; units of erg s _1 kpc -2 ) Sp a and S24 in 
the Pa recombination line and the MIPS 24 /im band, 
respectively, using a fixed photometric aperture matched 
to the MIPS 24 /jm resolution. The aperture size corre- 
sponded to a linear scale of ^300 pc at d = 5 Mpc. 

Among the high-metallicity H II knots, the range of 
LSDs observed by C+07 was 38 < log Sp^corr ^ 40 and 
40.5 < logS'24 < 42.5. An H II region with the same L24 
as M17, located in a galaxy 5 Mpc distant, would have an 
observed LSD of logS'24 = 41.00, well within the range. 
We can calculate the luminosity Lp a in the Pa line from 
the thermal radio continuum f iJ4.1.2[) . expressed as: 



Mff) 



(6) 



where j (Pa) and j,Aff) are the Pa and free-free emis- 
sivities |Osterbrocld 119741: IGilesI 119771 ) and F v {v) is the 
free- free radio continuu m flux density at frequency v 
(see IPovich et all [20071 for more detail). This gives 
Lp a = 4.1 x 10 37 erg s _1 . Assuming again d = 5 Mpc 
and a 300 pc test aperture, logSp Q = 38.76, also within 
the range of values reported by C+07 (since it is based on 
radio observations, there is no need to correct the LSD 
for extinction). 



Although M17 is a relatively low-luminosity H II re- 
gion in the extragalactic context, there is nothing un- 
usual about its LSDs compared to the C+07 sample. An 
M17 analog would be observable in the SINGS (Spitzer 
Infrared Nearby Galaxies Survey; iKennicutt et al.ll2003l) 
24 /im images as an H II knot in a nearby galaxy, assum- 
ing that it were not confused with a neighboring, brighter 
H II region, and it would fall well within the LSD distri- 
butions measured by C+07 for high-metallicity H II re- 
gions. 



4.2. Comparison of Mid-IR Galactic and Extragalactic 
SFR Calibrations 

We now generalize the analysis of M17 described above 
to a sample of eight Galactic H II regions (Tabled]). With 
one exception, this sample consists of H II regions that, 
like M17, meet all of the following criteria: (I) ionized by 
at least one early O-type star (spectral type 05 V or ear- 
lier) ; (II) sufficiently embedded that the bulk of the bolo- 
metric luminosity is reprocessed by dust and emitted in 
the IR; (III) ionizing cluster(s) with well-sampled (com- 
plete to m < 1.5 M©) XLFs reported in the literature; 
and (IV) evidence for ongoing star formation (for exam- 
ple, a detected population of stars with IR-excess emis- 
sion from disks). The exception is the ONC, which lacks 
early O-type stars but was included because it serves as 
an observational touchstone for massive cluster studies in 
general, and XLF studies in particular. Table [2] contains 
data relevant to SFR determinations for each H II region: 

Mass of stellar population. The total number of stars 
N pop associated with the ionizing cluster(s) was adopted 
from published XLFs (see references listed in Table [5]) 
and converted to total stellar mass, M pop , for m > 
0.1 M & using the Kroupa IMF. Uncertainties on M pop 
reflect only the statistical uncertainties on N pop , which 
include both the uncertain distance to each cluster and 
uncertainties introduced in extrapolating the XLFs to 
correct for incompleteness. Because we cannot always 
be certain that the XLF of a given cluster is complete, 
it may be safest to regard our reported Af pop values as 
lower limits. 



Star formation timescale and SFRsc- Following our 
selection criterion (IV) above, we assume that star for- 
mation has been continuous over the lifetime of each 
H II region, defined as the star formation timescale tsf. 
While X-ray emission gives no useful constraint on the 
ages of the stars, follow-up optical and IR photometry 
of the X-ray detected stellar population does. Values of 
tsf for each region are based on published results from 
isochrone fitting to optical/near- IR color-magnitude di- 
agrams (CMDs) and/or estimates of circumstcllar disk 
fractions using observed IR excess emission. Age deter- 
minations for very young (~1 Myr) stellar populations 
in regions with ongoing star formation are notoriously 
imprecise, since the observed age spreads are generally 
comparable to the measurement uncertainties and the 
reported ages themselves. We therefore use the upper 
end of the range of ages reported in the recent litera- 
ture (see references listed in Table [2J for each cluster. 
We had three motivations for this choice: (1) The diffi- 
culty in correcting for differential reddening in ncar-IR 
CMDs makes it easier to determine an upper age enve- 
lope than an average age for an obscured cluster. (2) 
Our values of Af pop may only be lower limits in some 



cases. (3) Our analysis of M17 ( tj4.1[) indicates that, in 
spite of a ^50% uncertainty on the age, SFRsc is signif- 
icantly higher than SFR based on diffuse tracers, hence 
we need only to place lower limits on the time-averaged 
SFRsc = Afpop/rsF (Table [2]) to test the generality of 
this result. We note that if there are bursts of higher 
SFR during the lifetime of an H II region, this would 
(temporarily) increase SFRsc- 

Mid-IR luminosity. Following iPovich et "all (|2007l) . 
we constructed IR SEDs for each H II region (Fig- 
ure Hh), using aperture photometry to extract 
background-subtracted fl ux densities from archival IRAS 
(jNeugebauer et abl 11984 ) and Midcourse Space Experi- 
ment (MSX; iPrice et al.ll200l survey images. We then 
performed a spline interpolation of the SED points, mea- 
sured the equivalent 24 /zm flux density F V (2A) as the 
mean value of the interpolation within the MIPS 24 /im 
bandpass, and calculated L24 using Equation [5j The 
dominant source of uncertainty in L24 is the uncertain 
distance to each region (Tabled]). 

Ionizing photon rates from radio continuum emission. 
We calculated N' c from the 5 GHz thermal radio contin- 
uu m flux of each H II region using the same conversion 
as iSmith et al.1 (|1978[ ) . Because these regions are large 
on the sky, most of the integrated flux measurements 
come from single-dish observations that are decades old. 
iPaladini et al.l (|2003l) compiled data from several large 
surveys of Galactic radio H II regions, and we use these 
cataloged flux densities for all regions except the Carina 
Nebula and the Rosette , for which w e adop t the values of 
ISmith fc Brooksl pOOl . and|Cdni]d (fl985l ). respectively. 
The available radio observations form a highly inhomoge- 
neous dataset, often suffering from background contam- 
ination and potentially unreliable calibration. We make 
no attempt to correct for internal absorption of ionizing 
photons and assume, as we found for the case of M17, 
that N c k, N' c . 

In Figure we plot logL 2 4 against logiV^, the re- 
lation that provides the foundation for the C+07 ex- 
tragalactic mid-IR calibration, for which N c is a direct 
proxy for SFR. We find a fairly tight correlation (rms 



scatter of 0.24 dex) with a marginally sub-linear slope 
(logi 2 4 oc [0.80 ± 0.06] log N' c ] dash-dotted line). The 
most significant outlier is the Rosette, which is the only 
"extended l ow-surface br ightness" radio H II region in 
the sample (Celnikl Il985[ ). hence it is possible that its 
measured radio continuum flux density includes a signif- 
icantly larger contribution from faint, extended emission 
compared to the other regions. C+07 derived a super- 
linear relation (log £24 cx [1.13±0.03] logA^; dashed line) 
for their extragalactic relation, and suggested that a lin- 
ear relation could apply to low-luminosity H II regions. 
Our data would be consistent with a linear relation if the 
ONC were excluded from the sample. 

In Figure [2]c we compare SFRsc with N' c . While H II 
regions with higher SFRsc also trend toward higher N{., 
all points (except for RCW 49 for which SFRsc is a lower 
limit only) lie above the linear correlation expected from 
the extragalactic calibration (Equation 01 plotted as a 
dashed line), some by as much as a factor of ~10 (dotted 
line). The weighted mean offset of the data (excluding 
RCW 49 and the ONC) from the predicted correlation is 
a factor of 2.2 (solid line), with very large scatter (0.38 
dex, rms). Note that the error bars on SFRsc reflect 
only the uncertainty in Af pop , which is not independent 
of the uncertainty in N' c (or L24 in Panel d) because 
both incorporate the distance uncertainty. The system- 
atic offset from Equation [4] might be explained in part by 
correcting N' c for absorption by dust and escape of ion- 
izing photons, but this correction is unlikely to shift any 
points by more than a factor of *~2, and as we showed 
for M17, some points would not shift at all. 

In Figure [Hi, SFRsc is plotted against L24; this corre- 
lation is significantly tighter than the above comparisons 
of SFRsc and N' c , with an rms scatter of 0.11 dex. The 
best-fit (dash-dotted; excludes RCW 49) line has slope 
0.76 ± 0.05, shallower than the 0.8850 slope of the anal- 
ogous C+07 extragalactic relation (dashed line). The 
C+07 calibration, scaled up by a factor of 2.7 (solid 
line), best matches the weighted mean correlation be- 
tween SFRsc and £24- Because our estimates of SFRsc 
(Table [2]) are lower limits, the factor of 2.7 represents 
a minimum systematic discrepancy between the mid-IR 
SFR calibrations and SFR measured from the resolved 
stellar populations. 

There are several possible explanations for the de- 
creased scatter in the SFRsc _ L24 correlation, as com- 
pared with the SFRsc — N' c correlation, including region- 
to-region variations in the correction factor for N' c to N c , 
inhomogeneity of radio observations, and/or stochastic 
sampling of the high-mass tail of the IMF. The last pos- 
sibility arises because ionizing photon rate produced by 
individual stars is a steeper function of mass than the 
bolometric luminosity, hence in extreme cases only the 
single most massive star in a given region dominates N c , 
while a wider range of stellar masses can heat dust and 
contribute to the 24 /mi emission. If the mass of the 
most massive star is vulnerable to large, stochastic fluc- 
tuations from cluster to cluster, we might expect L24 to 
trace SFR more smoothly and with less scatter than N c . 

In summary, we find that SFRsc derived from low- 
and intermediate-mass star counts to be systematically 
higher, by factors of at least 2-3, compared to either the 
calibration of Equation [4] (Figure [2]ci) or the C+07 extra- 
galactic SFR-L24 calibration (Figure [2]cQ . The discrep- 
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Fig. 2. — (a) Mid-IR SEDs (dots) for 8 Galactic H II regions, from archival IRAS and MSX data. Solid curves are spline interpolations to 
the SEDs. Effective luminosities in the Spitzer/MlPS 24 /im bandpass (L24; between the vertical dash-dotted lines) are plotted as boxes. 
(6) Plot of log L24, versus logA^, the uncorrected ionizing photon rate. The best fit to the data is plotted as a dash-dotted line with slope 
0.80 ± 0.06. C+07 derived a super-linear relation (dashed line; slope 1.13+/-0.03) for their extragalactic data, (c) SFRs derived from star 
counts, logSFRsCi plotted against logTV^. The relation given in Equation 4, a standard extragalactic calibration, is marked as a dashed 
line. Versions of Equation 4 scaled up by factors of 2.2 (fitting the weighted mean of the data, excluding RCW 49 and the ONC) and 10 
are plotted as solid and dotted lines, respectively, (d) Plot of logSFRsc versus £24- The correlation between logSFRsc and logZ/24 is 
consistent with the power-law relation derived by C+07 (dashed line) scaled up by a factor of 2.7 (solid line). The best-fit relation to our 
data is marked as a dash-dotted line and has a slope 0.76 ± 0.05. 



ancy between SFRsc and SFR24 (or equi valent ly SFR//) 
we derived for the specific case of Mf 7 fi M.1.3j) therefore 
appears to be a general feature of Galactic H II regions. 
Parallels might be drawn be t ween these results and the 
findings of iHeiderman et all (|2010() . who reported SFR 
surface densities measured from YSO counts in low-mass, 
Galactic molecular clouds that are higher by factors of 
5-10 compared to the standard e xtragalactic Schmidt 
Kennicutt Law (|Kennicuttlll998b[ ). 

4.3. Systematics Affecting SFR Calibrations Based on 
Lyman Continuum Photon Rates 

From the analysis in the preceding subsections, we 
conclude that even if Lyman continuum photon rates 
were known to arbitrarily high precision, an intrinsic, 
systematic uncertainty in the derived SFRs would re- 
main. Specifically, it appears that the mid-IR SFR cal- 
ibration of C+07 and, by extension, other calibrations 
that rely upon Lyman continuum photon rates underes- 
timate SFRs by at least a factor of 2-3 (Figures [3]c,d). 
This systematic affects the coefficient in Equation 01 the 
transformation between SFR and Lyman continuum pho- 
ton production rate. The coefficient can be parameter- 



ized as Mp 0p / (tN c ), where M pop /N c is the mass of a 
young stellar population normalized by its production 
rate of Lyman continuum photons, and the timescale r 
is either the duration of star formation traced by the 
observations of N c or the lifetime of the ionizing stars, 
whichever is shorter. M pop /N c is itself a function of 
age in p opulation synthesis models ( e.g., the Starburst99 
models; iVazquez fc Leithererl 120051 ) . but its functional 
form cannot be realistically modeled for cases where the 
SFR varies on timescales shorter than O-star lifetimes. 
As we cautioned in §2, Equation 4 strictly describes the 
SFR required to maintain a steady-state population of 
ionizing stars, implicitly assuming r « 5 Myr. Several 
effects likely to be important to M pop /N c and one impor- 
tant consequence of the assumed r are discussed below. 

4.3.1. Uncertain Lyman Continuum Photon Production 
Rates for O Stars 

Fundamentally, M pop /N c = (I/77) Y^,i( m /Qo)i, where 
m/Qo is the ratio of stellar mass to Lyman continuum 
photon production for individual O stars and ij is the 
(small) fraction of M pop contained in O stars according to 
the assumed IMF. Values for Q must be obtained from 
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models of O stars as functions of spectral type and lumi- 
nosity class . The Starburst99 grid of p opulation synthe- 
sis models (| Vazquez fc Leithererl 120051) incorporates the 
m odern O star calibr atio ns and synthet ic atm ospheres 
of lSmith et all (pool ) and iMartins etall (pool . These 
state-of-the-art calibrations employ non-LTE stellar dy- 
namic atmosphere models that include the effects of line- 
blanketed winds and predict cooler effective tempera- 
tures by as much as ^5000 K for the earliest spectral 
types. This results in ~40% reductions of Qn for a given 
spectral type com pared to earlier models (|Panagialll973t 
iVacca et al.l 119961 ) . Despite the obvious improvements 
in models of massiv e stars, we reiterate the caveat of 
IMartins et al.l (|2005l ): the uncertainty on Qo remains a 
factor of ~2 in any model, hence its impact cannot be 
neglected, although it is unclear in which direction this 
uncertainty might bias SFR determinations. 

4.3.2. The O star mass discrepancy 

The uncertainty in the numerator of m/Qo is also im- 
portant. Model-based calibrations of O star parameters 
predict the stellar mass as a function of spectral type in 
two different ways. The emergent spectrum gives values 
for stellar luminosity L, effective temperature T c g, and 
effective surface gravity g c g measured from spectral line 
broadening, giving 



Ms, 



(7) 



(jVacca et al.l 119961: IMartins et all 120051 ) where a is the 
Stcfan-Boltzmann constant and G is the gravitational 
constant. Alternatively, L and T c ff from stellar evolu- 
tionary tracks can be used to determine an evolutionary 
mass M BVO i, which is a ssociated with a stellar gravity 
gevoi (jVacca et al.l [l996). In general, g c s < g C voi, because 
a strong stellar wind lifts the surface layers of a star, re- 
ducing the gravity probed by spectral lines formed in the 
photosphere. Consequently, there is a discrepancy in the 
derived masses of up to a factor of ~2, a systematic effect 
with M spcc < M cvo \. O star masses based on M spoc in 
stellar models could therefore be undere stimates, par- 
ticula rly for the earliest spectral types ([Martins et al.1 
2005). This would tend to bias m/Qo downward in pop- 
ulation synthesis models, le ading to underes t imate s of 
M pop /N c and hence SFRs. iWeidner fc Vinkl (|20ToT) re- 
cently found good agreement between M spcc and dynami- 
cal masses calibrated to massive eclipsing binary systems, 
so it is possible that the O star mass discrepancy has been 
resolved. 

4.3.3. The Upper Mass Limit of the IMF 

The most massive stars (m > 40 M ) dominate N c 
but represent only 77 = 0.07 of A/ pop (assuming the 
Kroupa IMF) and follow, very approximately, the rela- 
tion Qo(m) oc m 3 5 . The ratio M pop /N c is therefore sen- 
sitive to the assumed upper mass lim it of the IMF. This 
cutoff is poorly known; for example, iZinnecker fc Yorkd 
(|2007h conclude that an upper limit of 1 50 M is con- 
sistent with data, while iCrowther et al.l (|2010l ) claim a 
maximum stellar mass of 300 M Q . This confusion may 
be compounded by mounting evidence that many O stars 
form i n binary systems wit h nearly equal-mass co mpo- 
nents ()Sana et al.l 120081 ). iMaiz Apellanizl ([2008;) per- 



formed numerical experiments testing the effect on ob- 
servational IMF determinations of mistaking unresolved 
massive binaries for single stars and concluded that while 
neglecting stellar multiplicity has a negligible effect on 
the high-mass IMF slope, the high-mass cutoff of the 
IMF will be systematically overestimated. In this work, 
we have conservatively assumed an upper mass limit of 
100 M (Equation [3]) . Increasing this limit would de- 
crease the SFR estimated from a given rate of ionizing 
photons, only worsening the discrepancy between SFR24 
and SFRso 

4.3.4. Runaway Massive Stars 

Studies of massive stars in the field (remote from any 
significant episode of star formation) show that ~10-30% 
of O stars escape their natal star-forming r egions (|Giesl 
fl!^IMoffat"enil[T998l: IQev fc Lambll20ill) . For young 
H II regions, O stars are likely ejected via dynamical en- 
counters, which are enhanced by their high binary frac- 
tion (e.g.. iLeonard fc Duncan] [19901: iKroupal l200l : di- 
rect signatures of O-star escape can be seen in imaging 
of bow shocks (jGva ramadz e et al.l 120111 ). Although the 
loss of ionizing stars from H II regions could contribute 
to the discrepancy between SFR24 and SFRgc , it can not 
account for the full magnitude of the discrepancy. 

4.3.5. The Assumed Form of the IMF 

Fundamentally, the discrepancy between SFRgc and 
SFR24 may indicate that M pop /N c predicted by stellar 
population synthesis models is systematically too low or, 
cquivalently, the fraction of mass 77 contained in mas- 
sive stars is too high. A straightforward way to increase 
M pop /N c would be to increase the numbers of low-mass 
stars (which contain the bulk of the stellar mass) rela- 
tive to massive stars (which dominate the luminosity in 
young populations). This is equivalent to steepening the 
intermediate- mass (or super-solar ) IMF slope. In their 
recent review, iBastian et al.l (|2010l ) argue that there is no 
convincing evidence for strong IMF variations, but this 
conclusion does not amount to a claim that the form of 
the "universal" IMF is known to high precision. As is 
the case with SFRs, IMF comparisons across different 
environments are more reli able than abs o lute m easure- 
ments. Indeed, Figure 2 of IBastian et al.l (|2010D reveals 
that while derivations of T in the super-solar mass regime 
are broadly consistent with the Salpeter-Kroupa slope 
(r = 1.3), this actually represents only a lower bound 
on r reported for 1.5 M < m < 5 M©. At interm ediate 
mass es, the IMF slopes cluster around T' — 1.7 (|Scalol 
[1991 . 

If we alter our standard Kroupa IMF (Equations 1-3) 
by increasing T = 1.3 to T' = 1.7 for 1.5 M < m < 5M , 
we derive SFR 24 = r)/r)'x SFR 24 = 1.4 x SFR 24 . We also 
find SFR^ C = 0.9 x SFR SC , by normalizing both IMFs to 
the same N pop . A steeper intermediate- mass IMF slope 
marginally reduces SFRs based on low/intermediate- 
mass star counts while significantly increasing SFRs de- 
rived from massive star diagnostics, driving the indepen- 
dent measurements toward convergence. The combined 
effect produces a factor of ~1.6 (— 1.4/0.9) increase in 
SFR24/SFRSC, enough to explain much of the systematic 
discrepancy. 

It would be reasonable to question the appropriateness 
of altering the IMF used in the XLF scaling method, 
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given that a fundamental assumption of this method is 
similarity with the ONC IMF ( U .1.11) Hence we note 
that the benchmark ONC IMF of lMuench et all ([20021 ) 
was measured using low-mass stars only (m < 2 M Q ; 
see Figure [1]), and T — 1.21 was simply extrapolated to 
the high-mass tail. In spite of its status as the near- 
est "massive" young cluster, the ONC contains too few 
intermediate-mass stars to support a statistical measure- 
ment of the super-solar IMF slope. We conclude that a 
super-solar IMF slope steeper than Salpeter is not ex- 
cluded by existing IMF determinations, and our method 
for comparing SFRs derived from XLF scaling and mas- 
sive star tracers can indirectly constrain its value. 

4.3.6. Stochastic Sampling of the IMF 

In the limit of very low SFR, the linear calibration 
of SFR against ionizing photon rate (Equation [4j must 
break down, as star formation can occur without pro- 
ducing any massive stars in a given region. In the case 
of most Galactic H II regions, we must consider the in- 
termediate possibility that the high-mass tail of the IMF 
will not be fully sampled. To representatively sample the 
high-mass tail of the I MF, a stellar population shoul d 
contain > 10 O stars (jCerviho fc Valls-OabaTil l2l)03Th 
which corresponds to Mp nn = 2, 9 00 M g (assuming a 
Kroup a or Chabrier IMF: lLee et al.ll2009f ) . TCervino et all 
(120031) also find that a total stellar mass of ~3, 000 M Q 
is the lowest mass for which sampling effects can be ig- 
nored (for a yo ung; stellar population and solar metallic- 
ity). Similarly. Ida Silva et al.l ([201 If ) find that, after ac- 
counting for stochastic effects, Lyman continuum fluxes 
are consistent on average with those expected for a fully 
sampled IMF for SFRs > 10~ 3 M Q yr" 1 (although the 
scatter in iV c /SFR is significant for SFR = 10~ 3 M Q 
yr -1 , approximately a factor of two). For SFR = 10~ 4 
Mq yr _1 , the average ionizing luminosity is suppressed 
by a factor of ~2. 

The eight Galactic H II regions described in Table [2] 
span a wide range of M pop and SFRs, with some sig- 
nificantly above these limits (Carina, M17) and others 
well below them (ONC, Rosette). In the lowest-mass re- 
gions, we might therefore expect calibrations based on 
population synthesis models with fully-sampled IMFs to 
underestimate SFR. However, stochastic sampling can- 
not explain the systematic offset of all H II regions in 
our sample, including the high-mass regions Carina and 
M17 that fully sample the IMF, from the C+07 relation. 
We note that the best-fit line to the SFRsc-^24 rela- 
tion (Figure [2ji) has a shallower slope than the C+07 
correlation, as would be expected if L24 preferentially 
underestimates SFRs in the lower-mass regions, such as 
the ONC. 

4.3.7. Star Formation Timescales 

As shown in Table [5J the (upper limit) ages for our 
Galactic H II region sample range from 1 to 5 Myr, with 
a median age of 2 Myr. These young ages are a conse- 
quence of our selection criterion (II), that most of the 
bolometric luminosity be reprocessed by warm dust and 
reradiated in the mid-IR. Older Galactic H II regions 
have been excluded from our analysis because they do not 
exhibit bright mid-IR nebulosity and would not be identi- 
fiable as 24 fim point sources in a SINGS image of an ex- 
ternal galaxy, hence they would not have been included in 



the C+07 calibration sample. Prominent excluded clus- 
ters include Westerlund 1, one of the most m assive young 
clusters in the Milky Way dClark et al . 2005) , and several 
Red Supergiant Clusters (|Figer et al.ll2006t iDavies et al.1 
I20071 lAlexander etail fcOOQ) , all of which contain signifi- 
cant populations of Wolf-Rayet stars, yellow hypergiants, 
red supergiants, and other evolved massive stars. The 
ages of these clusters fall in the 5-20 Myr range, hence 
they represent the future evolution of the giant H II re- 
gions in our sample. The dust in such evolved H II re- 
gions has either been destroyed by the hard radiation 
field or blown away by stellar winds and possibly super- 
novae, removing the source of the bright mid-IR emis- 
sion. Such regions certainly contribute to the extended 
IR emission of the Galaxy, since very little of their bolo- 
metric luminosity is reprocessed locally by dust. The in- 
tegrated Galactic IR emission is dominated by cool dust 
with a significant compon ent coming from older stars 
(e.g. iMisiriotis et al.l [20061. whereas the SEDs of young 
Galactic H II regions are dominated by warm dust heated 
by O stars (Figure [SJa) . 

The median value for the timescale tsf used to calcu- 
late SFRsc is a factor of ~2.5 times lower than O-star 
lifetimes, the timescale assumed implicitly by the C+07 
calibration and other calibrations based on population 
synthesis models, like Equation 0] (see caveat in §??). 
Hence the discrepancy between SFRsc and L24 (or N^,) 
in very young regions could be resolved by applying this 
correction factor to the calibrations of SFR versus diffuse 
emission tracers. 

The C+07 H II region sample was selected on the 
basis of compact, bright "knots" of mid-IR emission. 
Based on the discussion above, this selection criterion 
likely imposed an age constraint of tsf 5 Myr, which 
means that absolute SFRs derived by C+07 and sim- 
ilar studies of individual extragalactic star-forming re- 
gions could be systematically underestimated by factors 
of >8 Myr/5 Myr = 1.6, where 8 Myr is the steady- 
state timescale assumed implicitly by population syn- 
thesis (Section 2) and 5 Myr is the maximum age of a 
mid-IR bright H II region. We note that the H II region 
models employed by C+07 neglected H II region expan- 
sion, stellar winds, and supernovae, without which H II 
regions could remain mid-IR bright throughout the life- 
time of the ionizing stars. This discrepancy due to as- 
sumed timescales should not affect studies like K+09, 
where SFRs are measured as averages over entire galax- 
ies (although care is needed in the case of dwarf galaxies 
and other systems subject to short bursts of star forma- 
tion). In the case of normal spiral galaxies like the Milky 
Way, it is usually a good approximation that the galaxy- 
wide SFR has been constant over the last ^100 Myr, and 
the calibration of Equation [4] should hold. Indeed, K+09 
find that IR traces SFR more weakly in their data than in 
the C+07 sample, consistent with the longer timescales 
traced by galaxy-wide Ha and IR emission. 

5. ARE GLOBAL GALACTIC SFR ESTIMATES 
CONSISTENT WITH EXTRAGALACTIC DIAGNOSTICS? 

As we have seen, SFR estimates depend on a wide 
range of assumptions; incorrect input assumptions can 
lead to potentially large systematic errors in absolute 
SFR determinations. In addition, we must keep in mind 
that different assumptions may be appropriate in differ- 
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ent observational regimes, complicating attempts to com- 
pare SFRs. For example, if one set of assumptions were 
valid in external galaxies (which are often measured as 
wholes) while a different set of assumptions held in the 
Milky Way (where star formation is often studied piece- 
meal, and always in projection), these inputs could lead 
to a discrepancy between SFRs in the Milky Way and 
other galaxies. 

As reviewed in [J3l significant effort has been invested 
in measuring the Galactic SFR, but if we were external 
observers located Mpc from the Galaxy, would we still 
measure a Milky Way SFR of -2 M Q yr _1 ? The an- 
swer to this question is critical for comparing the Milky 
Way with other galaxies and extrapolating the detailed 
knowledge we gain in our Galaxy to more distant sys- 
tems. Here we present two tests that can assess, directly 
and with a minimum of assumptions, whether the Galac- 
tic SFR is consistent with extragalactic SFR calibrations. 
Although preliminary, these tests can serve as valuable 
sanity checks for our comparison of the Galactic SFR 
with calibrations used in other galaxies. 

5.1. The Brightest Radio Supernova Remnant 
( Cassiopeia A) 

The luminosity functions (LFs) of radio SNRs have 
been observed in ~20 galaxies and can be modeled as 
power laws with constan t index and scaling that de pends 
roughly linearly on SFR ([Chomiuk fe Wilc ots 2009| here- 
after CW09). Unfortunately, it is very difficult to es- 
tablis h a complet e sample of SNRs in the Milky Way 
le.g.. lGreeril27)n5h . so a comparison of the full LF is not 
a plausible technique for estimating the Galactic SFR. 
However, CW09 showed that there is a good correlation 
between a galaxy's SFR and the radio spectral luminos- 
ity of its most radio-bright SNR (L max ). We can use this 
£max~SFR correlation to make an independent estimate 
of the Galactic SFR, assuming that Cas A is the most 
luminous SNR in the Milky Way. This is a valid assump- 
tion because any SNR with a similar luminosity to Cas A 
would be bright enough for inclusion in any modern cat- 
alog, even if i t were located on the far side of the Galaxy 
(|Greenll200l . The details of our extragalactic sample 
are described in the Appendix; SFRs are measured us- 
ing the K+09 calibration of Ha and total IR emission. 
Here we measure SNR luminosities at 4.85 GHz because 
this frequency is less susceptible to the free-free absorp- 
tion that can dampen lower-frequency radio continuum 
emission in starburst galaxies. 

There is a clear correlation between L max and SFR 
(Figure E]), best fit at 4.85 GHz as: 

4i5 = (85if 7 ) SFR" 1±0 - 15 (8) 

This fit is shown as a solid line in Figure El This correla- 
tion is consistent with stochastic sampling of the power- 
law SNR LF (measured as described in the Appendix and 
CW09), marked as a dashed line in Figure El The scat- 
ter in the L max -SFR correlation is ±0.51 dex, consistent 
with the scatter expected from stochastic sampling of a 
power law (±0.55 dex). 

We note that IC 10 is a clear outlier (marked as 
a filled cir cle in Figure EH); its l arge "superbubble" 
(> 130 pc lYang fe Skillmanl [1993) is much more lu- 
minous than expected for a galaxy with IC 10's SFR. 
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Fig. 3. — For each galaxy in our sample, we plot the 4.85 GHz 
radio luminosity of its brightest SNR as a function of the galaxy's 
SFR. The solid line marks the best fit to the data. The dashed 
line represents the expected correlation due to random statistical 
sampling of a power-law LF, and the grey shaded region shows the 
expected la variation in the correlation due to stochasticity. The 
point marked by a filled circle is the outlier IC 10. The Milky Way 
is represented by a black star, assuming the Milky Way's most 
luminous SNR is Cas A and its SFR is 1.9 ± 0.4 M yr" 1 . 

lYang fe Skillm an (1993) claim that the most likely ex- 
planation for this source is multiple SNRs overlapping 
in the vicinity of a giant H II region. In fact, many 
of the brightest SNRs in galaxies may not be typi- 
cal of the SNR population as a whole, and they are 
likely to be a r ather heterogeneous class of objects 
(IChu et al.l fl999h . Some may be th e combined emis- 
sion from several overlapping SNRs ([Yang fe Skillmanl 
I1993D or the remnants of hyp ernovae/gamma-ray bursts 
( Lozinskava fe Mois cev 2007). Others may be interact- 
ing with unusually dense circumstellar material blown off 
by th e progenitor star (e.g., the l uminous SNR in NGC 
4449; Mili saylievic fe FeserJl2008D. Others may not be 
SNRs at all— IPakull fe Grisel (|2008l) point out that many 
large-diameter (200-500 pc) bright radio bubbles are 
likely blown by the jets of ultraluminous X-ray sources 
(ULXs) and not by SNe. For example, i n the most lu- 
mino u s radio "SNR" i n NG C 7793 (S26; IPannuti et all 
I2002H . IPakull fe Grisel (|2008D uncover a jet-like structure 
via high-resolution X-ray imaging which spans the ^450 
pc radio bubble and implies that that this bubble is 
blown by a microquasar and not by SN activity. Regard- 
less of the physical mechanisms driving the L max -SFR 
relation, it appears to hold empirically across more than 
three orders of magnitude, and we can use it as a test of 
current estimates of the Galactic SFR. 

In Figure El we plot the luminosity of Cas A and a 
Galactic SFR of 1.9 ± 0.4 M Q yr" 1 (as found in §3), 
marked with a solid star. The Milky Way appears to 
deviate from the L max -SFR relation as the most outly- 
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ing point besides IC 10. Using the technique described in 
the Appendix, we calculate a 0.2% probability that, given 
the luminosity of Cas A, the SFR of the Milky Way is 
1.9±0.4 M yr _1 (or, a 3% probability that it it > 1 M 
yr -1 ). Given the inherent uncertainty in the L max -SFR 
relation and the fact that outliers like IC 10 do exist, we 
can not conclusively eliminate the possibility that Galac- 
tic SFR estimates and extragalactic SFR diagnostics are 
consistent with one another. However, it is worth not- 
ing that the observed offset of the Milky Way from the 
L max -SFR relation implies that either the Galactic SFR 
has been overestimated or that extragalactic SFRs have 
been systematically underestimated. 

5.2. Luminous Young X-ray Point Sources 

We can further test if the measured Galactic SFR is 
consistent with SFR determinations in other galaxies us- 
ing young X-ray point sources (YXPs). The X-ray emis- 
sion from galaxies has often been claimed to trace SFR 
(eg .IGrimm et al.|[200l iRanalli et ai][200l IMineo et al.l 
120111 ). because of the high luminosities from high-mass 
X-ray binaries (HMXBs), in addition to SNRs and ULXs. 
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Fig. 4. — The integrated 2-10 keV luminosi ty of young X-ray 
point sources in each of 20 nearby galaxies {Persic & Rephaelil 



2007) plotted against SFR (assuming 50% uncertainty on Lyxp)- 
The best fit correlation is marked as a solid black line. The Milky 
Way is plotted as a black star, assumin g the SFR from §3 of 1.9 ±0.4 
M yr -1 and HMXB luminosity from lGrimnTet al. (2002). 



iGrimm et all (|2002| ) studied what the Milky Way X- 
ray source population would look like to an outside ob- 
server and estimated a total luminosity of Lhmxb = 
(2± 1) x 10 38 erg s -1 over 2-10 keV for HMXBs. We also 



add the luminosity of the Crab Nebula (Lx = 1 X 10 38 
erg s -1 assuming a dista nce of 2 kpc; IGrimm et al.ll2002l ; 
iDavidson fc Fesenl 119851 ). which dominates the Milky 
Way's SNR luminosity in the X-ray band. The Milky 
Way does not host any ULXs, so its Lyxp — (3 ± 1) x 
10 38 erg s _1 . 

We can compare the Milky Way with the sam- 
ple of 20 "local norma l and starburst" galaxies in 
iPersic &: Rephaelil (|2007f ). for which there are X-ray 
point-source LFs in the literature. Persic & Rephaeli ac- 
knowledge that the major uncertainty in using the X-ray 
point-source luminosity as a SFR indicator is contami- 
nation by low-mass X-ray binaries (LMXBs), which are 
associated with older stellar populations and trace stel- 
lar mass rather than SFR. In the case of the Milky Way, 
the integrated luminos ity of LMXBs is ten times greater 
than that of HMXBs; IGrimm et all 12002ft . Fortunately, 
the relative number of HMXBs as compared with LMXBs 
should b e a function of SFR d ivided by a galaxy ' s stel - 
lar mass (jGilfanov et al.l2004D . IPersic fc Rephaelil ((20071) 
deal with LMXB contamination by assuming that X-ray 
point-source LFs in systems with high SFR are repre- 
sentative of the young component and universally ap- 
plicable; similarly, they assume that LFs measured in 
E/S0 galaxies represent the old population of X-ray point 
sources. They fit each galaxy's LF with a combina- 
tion of old and young components. In Fig. [4j we plot 
the product of their measured point-source luminosi- 
tie s and YXP fract i ons, r escale d to the distances listed 
in iKennicutt et ail (|2008l ) and iMoustakas fc KennicuttJ 
(|2006D . We calculate SFRs as described in the Appendix, 
using the K+09 calibration. 

In Fig. HI we see that there is indeed a correlation be- 
tween Lyxp and SFR. The Milky Way, assuming a SFR 
of 1.9 ± 0.4 M yr -1 , is marked as a star. The stan- 
dard deviation around the best fit line is 0.47 dex, and 
the Milky Way is located ~1.5cr from the correlation. 
The Milky Way's deviation corresponds to an unusually 
high SFR for its measured Lyxp, an offset in the same 
directi on as observ e d for the L max -SFR relation. Re- 
cently, IMineo et al.l (|2011[ ) pointed out that the Persic & 
Rephaeli correlation is offset downwards from other X- 
ray binary-SFR relations in the literature by a factor of 
~ 2 — 2.5, fitting unusually low Lyxp values as a function 
of SFR. They hypothesize that IPersic fc Rephaelil (|2007l ) 
have overestimated the contribution from the old popu- 
lation of X-ray sources; if this is true, it would imply that 
the discrepancy of the Milky Way from the Lyxp~SFR 
relati on is larger t h an th at shown in Fig 2) In addi- 
tion, iLehmer et ail (|2010D show that the Milky Way is 
a slight (~ ler) outlier from their X-ray calibration to 
external galaxies. T hey find that given th e Milky Way's 
mass (5 x 10 10 M g : iHammer et al.ll2007h and SFR (as- 
sumed to be 2.0 Mq yr -1 , in good agreement with that 
used here), the ratio of X-ray luminosity in LMXBs as 
compared with HMXBs should be ~1.4, almost an order 
of magnitude lower than the Lt.mxb/Lhmxr ~ 1 mea- 
sured for the Milky Way bv IGrimm et al.l (|2002D . This 
is consistent with a SFR of 2.0 M Q yr -1 for the Milky 
Way being higher than expected from extragalactic cali- 
brations. 

5.3. Is the Milky Way an Outlier? 
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The Milky Way appears to deviate from both the 
-kyxp-SFR and L max -SFR relations at the 1 — 3a level, 
consistently showing a higher SFR than expected from 
the extragalactic relations. This deviation is marginal, 
and the Milky Way's similar offset from both correlations 
might simply be a coincidence, so we can not eliminate 
the possibility the the Milky Way's measured SFR is con- 
sistent with extragalactic determinations. However, both 
the SNR and YXP tests hint that a Galactic SFR of 2 
M.q yr _1 may be higher than we would measure for the 
Milky Way if we were external observers using the K+09 
SFR calibration. If this is true, there are three possible 
explanations for the difference: (1) estimates of SFR are 
biased high in the Milky Way (compared with "absolute" 
SFRs) (2) estimates of SFR are biased low in external 
galaxies; and/or (3) the Milky Way is a fundamental ly 
outlying system (as suggested bv lHammer 

If this deviation is real, it is unlikely to be related to 
the differences between Galactic H II regions and extra- 
galactic calibrations described in £14.21 Assuming that 
the most luminous SNR and the YXP population roughly 
trace the counts of massive stars, and understanding 
that the SFRs used in the L Y xp-SFR and L max -SFR 
relations are predominately measurements of the Lyman 
continuum rate, the Milky Way's displacement from the 
-^yxp-SFR relations implies that available mea- 
surements of N' c and L24 overestimate the Galactic SFR, 
whereas we have found that these quantities underesti- 
mate SFRs in individual Galactic H II regions. If the 
offsets described in §4.21 are due to a poorly constrained 
(but universal) IMF shape, they will affect the Milky 
Way and external galaxies equally, while differences due 
to assumptions about star formation timescales should 
be irrelevant when we average over entire galaxies (in- 
cluding the Milky Way). 

It is therefore difficult to explain the potential offset of 
the Milky Way from the L max /LYXP-SFR relations, but 
we must keep in mind that the Galactic SFR remains 
quite uncertain. For example, to calculate a total ioniz- 
ing luminosity for th e Milky Way from t h eir ca talog of 
WMAP H II regions, iMurrav fc Rahman! (|2010f ) had to 
make at least three significant assumptions, each of which 
increases the Galactic SFR by ^40-100% with large un- 
certainty: 1) they multiply their measured N' c by two to 
correct for distant H II regions missing from their catalog; 
2) they multiply by 1.5 to account for the diffuse compo- 
nent of the WMAP free-free emission, assuming this dif- 
fuse component has the same distribution as the detected 
H II regions; and 3) they multiply by 1.37 to correct for 
Lyman c ontinuum escape and abso rption by dust (detail 
given in iMcKee fc Williams! 119971 note that the K+09 
extragalactic calibration should also account for dust ab- 
sorption, because it uses a combination of Ha and IR 
emission) . In addition, if 1 /3 of the f ree- free emission has 
a very diffuse, smooth distribution as lMurrav fc Rahman! 
(2010) assert, some of this low-surface brightness emis- 
sion may go undetected in external galaxies, leading to an 
underestimate of extragalactic SFRs. Additional strate- 
gies for comparison between the Milky Way and external 
galaxies, along with better statistics on the imax/^YXP - 
SFR correlation, are required to test the provocative im- 
plication that the Galactic SFR may be overestimated, 
and/or extragalactic SFRs may be underestimated, by a 
factor of a few. 



6. CONCLUSIONS 

After normalizing published estimates of the Galactic 
SFR to Starburst99 population synthesis models and a 
Kroupa IMF, measurements converge to a Galactic SFR 
of 1.9 ±0.4 M yr" 1 . 

We test to see if SFR estimates in the Milky Way are 
consistent with extragalactic SFR diagnostics, using lu- 
minous SNRs and YXPs; the results are noisy, but the 
Milky Way is a 1 — 3 sigma outlier in both tests, implying 
that perhaps the SFR of the Milky Way is overestimated 
and/or extragalactic SFRs are systematically underesti- 
mated. 

We also test for consistency between SFR diagnos- 
tics within Galactic H II regions, where we can obtain 
exquisite detail and actual star counts. Our compari- 
son of SFRs based on tracers of massive stars against 
SFRs determined from star counts reveals that large sys- 
tematic errors persist in absolute SFR determinations. 
In particular, measurements of SFRs in both the Milky 
Way and external galaxies could be underestimating ab- 
solute SFRs by factors of at least 2-3. Two outstanding 
issues may explain the discrepancy between SFRsc and 
SFRs measured from indirect tracers: 

1. Uncertainty in the form of the intermediate-mass 
IMF. The systematic offset would be assuaged if 
the IMF has a slope that is steeper than Salpeter 
across the super-solar mass range (1.5 M Q < m < 
5 Mff l a tweak that is fully consistent with current 
data (jBastian et aTll2010H . 

2. Systematically overestimated star-formation 
timescales for extragalactic HII regions. The 
population synthesis models employed in extra- 
galactic SFR determinations often assume that 
the duration of star formation is longer than an 
O star lifetime. This assumption breaks down in 
the case of H II regions that are mid-IR bright, 
because such regions are systematically younger 
than the O-star lifetimes by factors of >1.6. 

We stress that these systematic issues are not mutually 
exclusive, and we cannot rule out additio nal contribu- 
tions from the O star m odels themselves (jSmith et al.l 
[20021: iMartins et alj[2005h and runaway O stars. In addi- 
tion, SFR determinations in the lowest-mass H II regions 
(< 3,000 M Q ) will also be affected by stochastic under- 
sampling of the high-mass end of the IMF. Because we 
derived only lower limits on SFRsc, the actual discrep- 
ancy may be higher than reported here. 

Nevertheless, distinguishing between the two primary 
scenarios is crucial to our understanding of absolute 
SFRs. In the IMF scenario, all Galactic and extra- 
galactic SFR determinations must be revised upwards, 
with far-reaching implications for applications that rely 
on absolute SFRs, like the Schmidt-Kennicutt relation, 
galaxy evolution models, and chemical enrichment by 
evolved low-mass stars. For example, if SFRs were 
larger on galaxy- wide scales, this would increase the 
discrepancy between SFR integrate d over cosmic time 
and measured stellar mass densities (Wi lkins et al]|2008t 
IChoi fc Nagarninil20ilh . Also ( assuming the Galactic 
SFR has been roughly constant throughout the history 
of the Milky Way), a larger amount of fresh un-processed 
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ISM material would be required to infall onto the Galac- 
tic disk to explain st ellar abundance pat terns and sustain 
star formation (e.g.. iSancisi et al1l2008l) . 

On the other hand, if the timescale scenario is suffi- 
cient to explain the discrepancy between SFRsc and the 
calibrations based on population synthesis models in a 
sample of IR-bright, H II regions, only SFR determina- 
tions in localized, very young star-forming regions will be 
affected, not entire galaxies. The extragalactic Ha+LR 
SFR calibration of K+09, for example, accounts for the 
ionizing luminosity of all massive stars in a galaxy, av- 
eraging across the star formation histories of many H II 
regions, so it should not be vulnerable to this systematic. 

Our discussion reminds us that implicit assumptions 
have the potential to significantly impact SFR determi- 
nations, both in a relative and absolute sense. In the fu- 
ture, additional tests and larger sample sizes are required 
to test the significance of the potential deviation of the 
Milky Way's SFR from extragalactic calibrations. Mean- 
while, additional multiwavelength observations of Galac- 
tic star-forming regions with a range of ages will have the 
power to discern between the two most likely scenarios 
for the H II-region SFR discrepancy. Such observations 
would allow for a more accurate determination of the 



IMF and also provide additional estimates of SFRsc for 
a range of star formation histories. This should pave the 
way to an independent, empirical calibration of absolute 
SFR versus tracers of ionizing luminosity. In order to 
place studies of galaxy evolution, ISM physics, and stel- 
lar populations on the most solid footing possible, we 
need to continue toward a better understanding of the 
systematics affecting SFR determinations, ensuring that 
Galactic and extragalactic calibrations achieve conver- 
gence. 
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APPENDIX 
DETAILS OF THE Lmax-SFR RELATION 

In an attempt to include as many galaxies as possible on the L max -SFR relation described in £ 15.11 we compile 
SNRs at 4.85 GHz, as opposed to the 1.45 GHz data used in CW09. Free-free and sychrotron-self absorption can 
greatly reduce the lower frequency (~l-2 GHz) radio continuum signal coming from S NRs in very dense , violently 
star- forming systems. However, frequencies around ~5 GHz remain unabsorbed (e.g., iParra et~atl [20071 ) . and are 
therefore ideal for searching for SNRs in starbursts. 4.85 GHz SNR samples are selected by the cri teria described in 
CW 09 and from the sam e references, with the addition of a few galaxies: NGC 1808 (radio data f romlSaikia et al.1ll990| 
andlCollison et al.lll994l ; excluding 58.56—36.3 which is proba bly the AGN nucleus of t he g alax y; Jimcn ez-Bailon et al.1 
2005|; NGC 4038/4039 (INeff fc Ulvestadl[2000l) : and Arp 299 (IPerez-Torres et al.H2009l and lUlvestadll2009D . In contrast 
with CW09, here we include the superbubble in IC 10. Although it is a severe outlier, being unusually luminous for a 
galaxy with IC 10's SFR, we cannot exclude it on these grounds alone as it might have an important impact on the 
imax-SFR relation. 

To calculate the luminosity of Cas A, we use a a 4.85 GHz flux d ensity of 806 Jy (jBaars et al.lll977l ). The luminosity 
of Cas A is fading with time (by ~1% per vea nlBaars et al1ll977l and these flux densities are standardized to epoch 
1980.0. We use a distance of 3.4 kpc to Cas A (|Reed et al.l 1995lT ^The fact that the luminosity of Cas A is decreasing 
implies that the Milky Way's position on the L max -SFR relation will change with time, but we note that this is also 
likely true for the other galaxies on the relation. This fact underscores the point that the L max -SFR relation is a 
statistical and stochastic relation; for it to hold true, as the most luminous SNR fades, another must be born and rise 
to take its place. 

Shortly after the analysis of CW09 was published, a new calibration for Ha+IR SFRs was released by K+09. Unlike 
the previous C+07 SFR calibration, the K+09 formulation is calibrated to entire galaxies, not individual star-forming 
regions. To calculate the extragalactic SFRs, we use the Ha+total IR (TIR) calibration from K+09, described as: 

SFR (M© yr" 1 ) = 5.5 x 10~ 42 (L Ha + 0.0024 L Tm ) (Al) 

where L# a and Ltir are both in units of erg s _1 , L# Q is corrected for foreground Galactic extinction, and Ltir is a 
linear combinat ion of luminosi t y mea surements in the Infrared Astronomical Satelli te (IRAS) 25, 60. and 1 00 /Ltm bands 
as suggested bv lDale fc Heloul (I2002D. Di s tance s and Ha measurements come from iKennicutt et alj (2008) if av ailable , 
and if not, fro mlMoustakas fc Kennicutti (|2006|). IRAS m e asure m ents come predomin antly f romlSanders et al.l ([2003), 
but also from [Fullmer fc Lonsdale! (|1989l ). iMoshir et al.l (|1992f ). iHunter et all (|1986D . and iRice et al.l (|1988l ). Basic 
information on the sample galaxies and derived SFRs can be found in Table [3] 

Table [3] also lists the number of SNRs (Nsnr) in the complete samples at 4.85 GHz. Many SNR surveys do not 
cover the entire extent of the galaxy, so we also include the field of view of the SNR survey (FoV), and the fraction of 
the SFR covered in this field of view (SFR Frac). The SFRs used in the L max -SFR relation are the measured SFRs 
times these fractions. 

As in CW09, we find that the SNR LF can be described as a power law: n(L) = A L@ . Using the techniques 
described in CW09, we find, for the 4.85 GHz data, a best-fit power-law index of j3 — —2.02 and near-linear scaling of 
the LF with SFR as: 

A 4 .85 = (75t\l) SFR 119±o n . (A2) 

These measurements of the SNR LF allow us to quantify the fact that galaxies with higher SFRs are more likely to 
host very luminous SNRs because they have larger populations of SNRs. 
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TABLE 3 

Galaxies Used in the L max -SFR Relation for SNRs 



Galaxy 


R.A. (2000) a 


Dec (2000) a 
(":':") 


Type a 


Distance 


SFR 


SFR Frac 


FoV 


NsNR 




(hr:min:sec) 




(Mpc) 


(M Q yr" 1 ) 




(arcmin) 


(4.85 GH: 


LMC 


05:23:34.5 


-69:45:22 


SB(s)m 


0.05 ± 0.003 


0.21 


1.0 


600 


26 


SMC 


00:52:44.8 


-72:49:43 


SB(s)m pec 


0.06 ± 0.003 


0.03 


1.0 


270 


11 


IC 10 


00:20:17.3 


+59:18:14 


IBm b 


0.66 + 0.03 


0.02 


0.91 


4.3 


4 


M 33 


01:33:50.9 


+30:39:36 


SA(s)cd 


0.84 ± 0.04 


0.25 


0.98 


10 


35 


NGC 1569 


04:30:49.0 


+64:50:53 


IBm 


1.9 + 0.4 


0.24 


1.0 


2.2 


18 


NGC 300 


00:54:53.5 


-37:41:04 


SA(s)d 


2.00 + 0.10 


0.10 


0.65 


9 


8 


NGC 4214 


12:15:39.2 


+36:19:37 


IAB(s)m 


2.92 + 0.29 


0.15 


0.77 


2.2 


3 


NGC 2366 


07:28:54.6 


+69:12:57 


IB(s)m 


3.19 + 0.32 


0.08 


1.0 


4.6 


4 


M 82 


09:55:52.7 


+69:40:46 


10 


3.53 + 0.35 


3.51 


1.0 


0.5 


23 


M 81 


09:55:33.2 


+69:03:55 


SA(s)ab 


3.63 + 0.18 


0.52 


0.88 


12.5 


3 


NGC 7793 


23:57:49.8 


-32:35:28 


SA(s)d 


3.91 + 0.39 


0.28 


1.0 


9 


4 


NGC 253 


00:47:33.1 


-25:17:18 


SAB(s)c 


3.94 ± 0.39 


2.89 


0.69 


2.75 


18 


NGC 4449 


12:28:11.9 


+44:05:40 


IBm 


4.21 + 0.42 


0.41 


1.0 


5 


5 


M 83 


13:37:00.9 


-29:51:56 


SAB(s)c 


4.47 + 0.22 


2.08 


0.98 


9 


14 


NGC 4736 


12:50:53.0 


+41:07:14 


(R)SA(r)ab 


4.66 + 0.47 


0.39 


0.81 


2 


10 


NGC 6946 


20:34:52.3 


+60:09:14 


SAB(rs)cd 


5.9 ± 1.2 


2.61 


0.99 


9 


21 


NGC 4258 


12:18:57.5 


+47:18:14 


SAB(s)bc 


7.98 ± 0.40 


1.33 


0.79 


7 


4 


M 51 


13:29:55.7 


+47:13:53 


SAbc 


8.00 + 1.04 


2.30 


1.0 


9 


30 


NGC 1808 


05:07:42.3 


-37:30:47 


(R)SAB(s)b 


10.6 + 2.1 


2.52 


0.60 


0.3 


10 


NGC 4038/39 


12:01:53.3 


-18:52:37 


Pec 


25.0 + 5.0 


10.4 


1.0 


2 


35 


Arp 299 


11:28:32.2 


+58:33:44 


Pec 


50.6 ± 10.1 


56.0 


0.5 


0.01 


19 


Arp 220 


15:34:57.1 


+23:30:11 


Pec 


85.2 ± 17.0 


77.2 


1.0 


0.02 


12 



a From NED. 6 From lde"Vaucouleurs fc Freemanl i 19721 ). 

We can quantify the probability that the Galactic SFR converged upon in is consistent with the by L max -SFR 
relation by assuming that: (a) the L max -SFR relation is the result of stochastic sampling of a power law; and (b) 
the power law scaling can be described by Equation IA2I For a given SFR, we then calculate how many SNRs can 
be expected for a power law with this scaling, and Monte-Carlo sample that number of SNRs from the power law 
distribution 10 4 times. This then provides 10 4 values of L max for the given SFR, and we repeat this process for a 
range of SFRs from 0.005 to 20 M Q yr -1 . For each SFR, we measure what fraction of these L max values fall within a 
factor of 1.5 of the luminosity of Cas A. This probability as a function of SFR has a roughly log- normal form, peaking 
at 0.2 M yr" 1 . 

If we use the 1.4 GHz data presented in CW09 (which has a smaller sample size than the 4.85 GHz data) for an 
analysis of the Milky Way on the L max -SFR relation, we find that the Milky Way is offset from the correlation in the 
same direction, although it is less aberrant (giving a 15% chance that the SFR of the Milky Way is 1.9±0.4 M Q yr -1 ). 



